Scalar Field Dark Energy Perturbations and their Scale Dependence 
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We estimate the amplitude of perturbation in dark energy at different length scales for a 
quintessence model with an exponential potential. It is shown that on length scales much smaller 
than hubble radius, perturbation in dark energy is negligible in comparison to that in in dark 
matter. However, on scales comparable to the hubble radius (Ap > lOOOMpc) the perturbation in 
dark energy in general cannot be neglected. As compared to the ACDM model, large scale matter 
power spectrum is suppressed in a generic quintessence dark energy model. We show that on scales 
Xp < lOOOMpc, this suppression is primarily due to different background evolution compared to 
ACDM model. However, on much larger scales perturbation in dark energy can effect matter power 
spectrum significantly. Hence this analysis can act as a discriminator between ACDM model and 
other generic dark energy models with Wde 7^ — 1- 

PACS numbers; 98.80.-k, 95.36.+X, 98.65.-r 



I. INTRODUCTION 



Cosmological observations suggest that about 70% 
content of our universe is made of a form of matter which 
drives the accelerated expansion of the universe [H. These 
observations include Supernova type la observations [2|, 
observations of Cosmic Microwave Background [1, 0,13 
and large scale structure 0, H, Q ■ The accelerated ex- 
pansion of the universe can of course be explained by 
introducing a cosmological constant A in the Einstein's 
equation [l^, However, the cosmological constant 

model is plagued by the fine turning problem p^ . This 
has motivated the study of dark energy models to ex- 
plain the current accelerated expansion of the universe. 
The simplest model as an alternative to cosmological con- 
stant model is to assume that this accelerated expan- 
sion is driven by a canonical scalar field with a poten- 
tial ¥{(/)). This class of dark energy models are known 
as quintessence models and the scalar field is known as 
a quintessence field. Various quintessence models have 
been studied in literature [H, [3 13 ■ There exists an- 
other class of string theory inspired scalar field dark en- 
ergy models known as tachyon models [l^ [l3| ■ Models 
of dark energy which allow w < —1 are known as phan- 
tom models Phantom type dark energy can also be 
realized in a scalar tensor theory of gravitation. (See for 
example Ref. 
essence field 



1£|.) Other scalar field models include k- 
15 1, branes [l^], Chaplygin gas model and 



its generalizations [2l| . There are also some phenomeno- 
logical models [l^, field theoretical and reorganization 
group based models (see e.g. \23l ). models that unify 
dark matter and dark energy |24l |. holographic dark en- 
ergy models [1^ and many others like those based on 
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horizon thermodynamics (e.g. see [2a|). For reviews of 
dark energy models see for instance Ref. [27*1 and for con- 
straining parameters using observations see Ref. 28] . 

Homogeneous dark energy distribution leads to accel- 
erated expansion of the Universe which, in turn, governs 
the luminosity distance and angular diameter distance. 
The rate of expansion also influences growth of density 
perturbations in the universe. This is evident from the 
abundance of rich clusters of galaxies and their evolution 
and the Integrated Sachs Wolfe effect [29j . 

In this paper we present a set of argument which lead 
to the conclusion that inhomogeneous dark matter with 
homogeneous dark energy at all length scales is inconsis- 
tent with Einstein's equation if pde 7^ —pde- We further 
analyze how dark matter power spectrum is influenced by 
perturbation in dark energy with an evolving equation of 
state parameter. 

Dark energy perturbations have been extensively stud- 
ied in the linear approximation ^33, Issl [3^ . It was shown 
in Ref. [33] that dark energy perturbations affect the low I 
quadrapole in the CMB angular power spectrum through 
the ISW effect. This analysis was done for a constant 
equation of state parameter. For models with w > —1 
this effect is enhanced while for phantom like models it 
is suppressed. In these models dark matter perturba- 
tions and dark energy perturbations are anti-correlated 
for large effective sound s pee ds. This anti-correlation is 
a gauge dependent effect [331 ■ Detailed studies in dark 
energy perturbations also include [H, [s^ . Clustering of 
dark energy within matter over density and voids were 
studied by Mota et al[31]. 

In this paper we study the evolution of perturbation 
in dark energy in a quintessence model which results in 
an evolving equation of state parameter different from 
that considered in Refs.fs^, HI]- We use a specific model 
of scalar field dark energy with an exponential poten- 
tial. We find that although for scales much smaller than 
hubble radius the perturbation in dark energy is small, 
for scales > the perturbation in dark energy can 
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be comparable to that in matter. Hence, although on 
small scales (<< H~^) the dark energy can be treated 
as homogeneous, one has to take into account the per- 
turbation in dark energy over scales ~ if Wde 7^ — 1- 
Clearly, in the specific case of Wde — the dark energy 
is homogeneous at all scales. 

We choose to work in the longitudinal gauge as in that 
case we can directly relate the metric perturbation $ to 
the gravitational potential perturbation. For a specific 
model, we investigate how quintessence dark energy in- 
fluences matter power spectrum. We show that on scales 
Xp < lOOOMpc, the matter power spectrum is not signifi- 
cantly affected whether or not we include fluctuations in 
the quintessence field in perturbation equations. How- 
ever, on much larger scales, including or excluding fluc- 
tuations in quintessence fleld does result in significant 
changes in the matter power spectrum. 

This paper is organized as follows. In Section |TT] we dis- 
cuss the background cosmology for matter and the scalar 
field system and describe the cosmological perturbation 
equation in longitudinal gauge for this system. In Section 
mil we obtain the numerical solution of the perturbation 
equation to determine the ratio of the perturbations in 
dark energy to the perturbations in matter at various 
length scales. Section ITVl summarizes the results. 



II. INHOMOGENEOUS MATTER AND DARK 
ENERGY 

We shall consider a system of minimally coupled mat- 
ter (dark -I- baryonic) and canonical scalar field ^ with 
Lagrangian of the form : 



(1) 



As a simple example, we consider a scalar field poten- 
tial of the form El, IH, IH : 



(2) 



Here A and Vo are two parameters of the potential and 
Mp = (87rG)-i/2 is the Planck mass. This potential leads 
to scaling solutions [13, 0, [3§| • However, for treating ex- 
ponential potential as a possible candidate for dark en- 
ergy, we require that the scalar field should not enter the 
scaling regime. This is possible to achieve by restricting 
the choice of the parameter in the exponential potential 

A spatially flat homogeneous and isotropic line element 
is described by the FRW metric of the form : 



(3) 



For the system of pressureless matter and quintessence 
with exponential potential the background evolution a{t) 
is completely determined by solving the following Fried- 
mann equation and the Klein Gordon equation : 



Pm.a 



and 



3^0-^l^„exp[-VAA].O 



(4) 



(5) 



A. Equations for Perturbations 

For analyzing perturbations in scalar fields and mat- 
ter we work in the longitudinal gauge (ill . IH, (For a 
recent pedagogical review see, l44|)- Fo'" scalar field and 
for pressureless matter (described as perfect fluid), the 
perturbed energy momentum tensor has no anisotropic 
term i.e. 6T^j is diagonal (where i,j — 1,2,3)^. Ein- 
stein's equations would then imply that the scalar metric 
perturbations in this gauge are completely described by 
a single scalar variable $. For such a system, perturbed 
FRW metric in longitudinal gauge attains the form [isj : 



(1 -t- 2$)dt2 - a^it){l - 2$)[da;2 -|- dy^ + dz^] (6) 



It is a good approximation to treat Dark Matter, bary- 
onic matter, etc. as perfect fluids. The energy momen- 
tum tensor of a perfect fluid is described as: 



= [p + Ph^^u, - pS^ 



(7) 



Perturbations in the energy density p, pressure p and 
the four velocity u'^ are defined as: 



p{t,x) 
p{t,x) 



Po{t) 
Po{t) 



5p{t, x) 
6p{t, x) 



(8) 
(9) 

(10) 



where = {1,0,0,0}; po{t) and po{t) are average 

values of the energy density and pressure respectively on 
a constant time hypersurface. The peculiar velocity is 
given by Su\ Substituting Eqs. dH]), ® and ^ in Eq. 
([7]) and neglecting second and higher order perturbation 
terms we get: 



dT% - Sp 



{po+Po)Su'- 

-Sp5^^ 



(11) 
(12) 

(13) 



^ We shall denote scalar field by {j> and the metric perturbation in 
the longitudinal gauge by $ 



^ Please note that through out this paper = 0,1,2,3 and i, j = 
1,2,3 
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For scalar (quintessence) field with Lagrangian of the 
form Eq.(IT]), we define the perturbations as : 



(14) 



where 0o(i) is the average value of the scalar field on 
the constant time hypersurface. The energy momentum 
tensor for the scaler field is given by 



(15) 



Substituting Eg. p^ in Eq. p3|) and subtracting the ho- 
mogeneous part in the energy momentum tensor we get 



B. Inhomogeneous Dark energy 



(16) 
5] (17) 
(18) 



The existence of inhomogeneity in non relativistic mat- 
ter (dark -f baryonic) is evident from direct observations. 
We ask : is it reasonable to assume that dark energy is 
in general homogeneous given the observational fact that 
matter (dark -I- baryonic) is clustered? Dark energy has 
equation of state p < — leading to the accelerated 
expansion of the universe. The fluid with such an equa- 
tion of state behaves gravitationally as a repulsive form 
of matter and opposes gravitational clustering. 

Our assumption that matter and scalar field are mini- 
mally coupled implies that energy momentum tensor for 
both matter and scalar field are individually conserved. 
This would then imply that for both these components 
the variables defining the perturbations Sp, Sp and Su' 
would satisfy the following two equations^ : 
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FIG. 1: This plot shows the variation of Qm and fi^ as a 
function of scale factor. We have chosen A = 1. In the x-axis, 
"Log" refers to the logarithm to base 10. 



an observational evidence that the matter in the universe 
is clearly not homogeneously distributed. Hence, if we as- 
sume that dark energy with ui ^ — 1 is homogeneously 
distributed at all length scales, then it is inconsistent 
with the observed features of the Universe fs^, [H, [13. 



C. Linearized Einstein's Equation 

The perturbed Einstein's Equation about a flat FRW 
metric are given by dG^ ^ = SttG 5T^^^,. In our case 
5T^, = 6T^^^^^,,^^^ + ST'^^^^y Since matter has neg- 

ligible pressure we set ST^-^^^^^^^-^ — 0. Fluctuation in 
pressure is contributed only by the scalar field. Calculat- 
ing the perturbed Einstein's tensor ^ from the line 
element dH) and substituting 5T^'^ from Eqs.(Ill|) to ^ 
and from Eqs. (fT6|) to (fT7|) . we obtain the following lin- 
earized Einstein's equations : 



5p+[po+Po)V.5u + ?,-{5p + 5p)^-i<^{po+Po) = (19) 
a 



(Po +Po)Su.i +Po6Ui - 6p^i - <i>,,;(po +Po) (20) 

If dark energy with w ^ —1 were to be homogeneously 
distributed, i.e if ST^^^^je) = 0, then, Eqs.dH]) and ^ 
imply that the gravitational potential does not depend 
on space and time (i.e. $ = constant). If this is the case 
then in the line element © , we can rescale time and space 
coordinate such that Eq.® becomes FRW metric 
This would then imply that 5pm = 0. Hence homoge- 
neously distributed dark energy at all length scales would 
naturally imply that matter is also distributed homoge- 
neously. As we see structures over different scales, this is 



3 Note that for scalar field ST" - = 
written as ^ = {pif,^ +P<j)o)<5'"i 



= (poSip^- This could also be 
(<,) where Su^^^^ = 4>~'^5(p^i . 



$ + 4-$+ (2- + ^)$ 



$ + 

a 



-AttG 



5pm + (j^oH- 



+ r (0o)<50 (21) 

47rG 

-^4>l-V'{<i>o)5<p\ (22) 



where V'{4)o) = dV{(j)o)/d(j)o and is the potential for 
the matter peculiar velocity such that Sui = ViVm- Since 
these equations are linear, we have fourier decomposed 
the perturbed quantities such as $, d(j), Spm and Vm and 
replaced by — fc^, where k is the wave number defined 
as fc = 27r/Ap and Ap is the comoving length scale of per- 
turbation. In these equations all the perturbed quantities 
correspond to the amplitude of perturbations in the /c*'' 
mode. 
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FIG. 2: This plot shows the variation of equation of state 
parameter as a function of scale factor. We have chosen 
A = 1. In the x-axis, "Log" refers to the logarithm to base 
10. 



FIG. 3: This plot shows the how the gravitational poten- 
tial evolve with time. The bold line shows how gravitational 
potential $jv(a) evolves with the scale factor in the AC DM 
model, while the dotted lines shows (a) for the dark energy 
model considered in this paper. Here the parameter A = 1. 
In this plot the wave length of perturbation Ap — 10^ AIpc. 



Eq. ([22|) is the dynamical equation for the metric per- 
turbation $ and the perturbation in the scalar field turns 
out to be the driving term. The unknown variables are $, 
6(j) and Spm- Once and S(j)(t) are known then the po- 
tential for matter peculiar velocity Vm can be determined 
using Eq. l|^ . Also it is interesting to note that once ^{t) 
and S(l)(t) are known then even the matter density per- 
turbation Spm can be determined using Eq. (|2ip . Hence 
for such a system of pressureless matter and scalar field, 
the dynamics of perturbations is uniquely determined if 
we know the solution and 6(j){t). Hence we need just 
two second order equations connecting 5(j){t). We 

choose Eq. ([^^ as one of these equations. In addition 
to this, the dynamical equation for the perturbations in 
the scalar field 6(p{t) is obtained from the scalar field la- 
grangian ([T]) and this is given by : 



S(P + 3-1 
a 



2$l/'(0o)-4$0o 



V"{<j)o)S^ - 
(24) 

For any quintessence potential V{(f>), in a system 
of matter and scalar field, one can solve the coupled 
Eqs.([22| and ((24|) to study the behavior of the perturbed 
system. Once the solution ^{t) and S(j){t) is obtained, 
we can then calculate the fractional density perturbation 
defined as: 



6 



Sp 



(25) 



for both matter as well as scalar field from Eqs. (fT6|) and 
(PT|) . This is given by : 



ll + V'iMSc^ (26) 



+ ^(</'o) 



^ 3^$ + 3-$ + \ 



Pra^a 



\4>l + v{<t>o) 



(27) 



Using these two equations we shall calculate 5„i and 
Srj, for a system consisting of dark matter with negligible 
pressure and quintessence dark energy with exponential 
potential. 



III. NUMERICAL SOLUTIONS 

For solving the background equations [Eqs.(|4]) and ([5])], 
we introduce the following dimensionless variables: 



y 



Mr> 



a 
a,: 



X = Hi{t-ti) 



(28) 

(29) 
(30) 



where a^, <j)oi and Hi are the values of the scale factor, 
scalar field, and the bubble parameter at some initial time 
t = U. 

In terms of these new variables, the two equations 
[Eqs. (HI and ([5])] describing the background cosmology 
becomes: 



y'2 



y 



s 



yexp(-VAy)j = (31) 
VAFexp(-%/Ay) = (32) 



where prime " ' " corresponds to the derivative with re- 
spect to X, = VoiJj^^M"^ exp(— i/AMp'^^oi) and 
is the dimensionless matter density parameter at the ini- 
tial epoch t — ti. 

Since f2totai = f^mi + f^<Ai = 1, it follows that: 



V 

y[ 



1 



(1 - VLmi){l - Wi) 



v/3(l-r!„„)(l + m0 



(33) 
(34) 
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FIG. 4: This plot shows how matter perturbation Sm evolves 
with scale factor. The bold lines on this plot corresponds 
to ACDM model, while the dotted lines in the same plot 
shows how the matter perturbations evolves with scale factor 
for the dark energy model considered in this paper. Here 
the parameter A = 1 and the wave length of perturbation 
Ap = 10^ Mpc. 



FIG. 5: The ratio of power spectrums P{k)QCDM /P{k)ACDM 
as a function of length scale of perturbation Ap. If the poten- 
tial parameter A = 1, then it turns out that at the present 
epoch, the equation of state parameter wo = —0.83, and if 
A = 0.1 then wq = —0.98. In the x-axis, "Log" refers to the 
logarithm to base 10. 



where Wi is the value of the equation of state parameter 
of the scalar field (j) at t = ti. Choosing Wi to be very 
close to —1 at say red shift of Zi = 1000, we solve the two 
background equations [Eqs.Q and dS])]. 

Fig[T] shows the plot of fi^ and Jl^ as a function of 
scale factor. This figure gives the value of scale factor 
at the matter dark energy equality aeq — 0.68. The 
corresponding redshift of matter dark energy equality is 
Zeg = 0.46. The red shift at which the universe under- 
went a transition from decelerated expansion phase to 
accelerated expansion phase turns out to be Zacc = 0.81. 
Figl2] shows the evolution of the equation of state param- 
eter W0 as function of scale factor. For the choice A = 1, 
this gives the value of the at the present epoch to be 
w^^ = -0.83. If A = 0.1, then w^^ = -0.98. 

For numerically solving the two perturbation equations 
Eqs.([m and ([M]) . we introduce the following two dimen- 
sionless variables : 



5y 



(35) 



(36) 



Here is the normalized gravitational potential with 
<l>i being the value of the metric potential at the initial 
time t = ti. In terms of these two new variables, the Eqs. 
([22|) and (|24|) can be rewritten as 



5y" + 3-<5y' + {\5y - VA)t/ exp(-VAy) - 

4'J>;y' + |J=0 (38) 

Here k = k/Ho, where Ho is the bubble parameter at 
the present epoch, is the wave number scaled with respect 
to the bubble radius. In Eq. ([55)1 the constant f3 = -pj, 

where sj, is the value of s' at the present epoch. This can 
be determined numerically from the Eq. (|3ip . 

We assume that the perturbation in the scalar field in 
the matter dominated epoch at z w 1000 is negligibly 
small compared to other perturbed quantities such as $, 
6m etc. The scalar field can then be treated as initially 
homogeneous at t = ti. This corresponds to setting the 
initial condition Syi ~ and Sy'^ = 0. The only initial 
condition that needs to be determined is the value of 
a.t t ~ ti. In the matter dominated epoch (z « 1000), 
the linearized Einstein's equation (|22p can be analyti- 
cally solved to obtain the solution <i>(t) cx t^^^^. Hence 
<i>(t) decays to zero in the matter dominated epoch for 
all values of wave number k and we can set the initial 
condition (f>^ (fc) = 0. 

After solving the two perturbation equations using the 
above initial conditions we can now determine the di- 
mcnsionlcss density perturbations Sm and defined in 
Eq. (P7|) and Eq. (pS)) . In terms of the dimensionless vari- 
able defined in Eq. ([55)1 and Eq. we can express Sm 
and Sijj in the following form: 



^ y'Sy' - <i>„y'^ - VXSyV exp{-VXy} 
" ^ + Vexp{-y/Xy} 



(39) 



2- + — 



y'Sy' 



+ VXV exp{-VXy)Sy] = 



(37) 



Sm 



S'2 



s' P 
2-$'.. +2-^$„ -I- 



{y'Sy' - - VXSyV eM-^y})] (40) 
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FIG. 6: This plot shows how perturbations in dark energy 5 4, 
evolve with scale factor. Here the parameter A = 1 and the 
wave length of perturbation Ap = 10^ Mpc. 
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FIG. 7: This plot shows the dependence of the ratio b^l&m at 
the present epoch on the wave length of perturbations Ap. In 
this plot the parameter A = 1. In both the axis, "Log" refers 
to the logarithm to base 10. 



In Fig. [3]we have plotted as a function of scale fac- 
tor for the value of the parameter A = 1. The dotted line 
in this figure shows ^ ^ (a) for the model of dark energy 
considered in this paper. We can see that in the matter 
dominated epoch = 1 and it is constant. Once the 
dark energy dominated epoch begins, the gravitational 
potential starts to decay. For comparison, in the same 
figure we had plotted [see bold line in Fig[3] the form of 
$„(a) in the ACDM model. This impHes that the be- 
havior of <i>„ in the matter dominated era is the same in 
both ACDM model and the model of dark energy consid- 
ered in this paper. But once the dark energy dominated 
phase begins, $„(a) decays faster than the correspond- 
ing (a) in the ACDM model. In this figure the wave- 
length of perturbation was fixed to Ap = lO^Mpc. With 
Ho — TiKm/s X Mpc~^\&], this value of Ap corresponds 
to fc = k/Ho ~ 26. 

In Fig. m we have plotted the matter perturbations 5m 
as a function of scale factor for the dark energy potential 
parameter A = 1. In this plot the solid line corresponds 
to (5m (a) for ACDM model and the dotted line in the 
same plot shows 5m{a) for the dark energy model consid- 
ered in this paper. Matter perturbation 5„i{a) initially 
grows linearly with scale factor in the matter dominated 
epoch, however, once the universe undergoes transition 
from the decelerated expansion phase to accelerated ex- 
pansion phase, the growth of 5m{a) is suppressed. In 
fact the growth of 5m (a) is suppressed more in the dark 
energy model with exponential potential than the cor- 
responding 5m{a) in the ACDM model. At the present 
epoch, for the same set of initial conditions we find that: 

^^•^^^(J^^Cz = 0) = 0.96 (^^^*^)<5^(z = 0) (41) 

where QCDM corresponds to the quintessence -f cold 
dark matter. Hence, the matter power spectrum is sup- 
pressed in QCDM model considered in this paper than 
the corresponding value in the ACDM model. This 
nearly 4% suppression corresponds to the length scale 
of perturbation Ap — \Q^Mpc and for the value of 
the parameter A = 1. By 4% suppression, we mean 




FIG. 8: This plot shows the dependence of the ratio 5q,/5m at 
the present epoch on the equation of state parameter Wo- It 
is the value of the parameter A in the potential V{(f)) which 
determine the present value of the equation of state parameter 
Wo- In this plot the wave length of perturbation Ap = 10'^ Mpc. 



that [(P(fc) ACDM — P{k)QCDM)/ P{k)ACDM] X 100 = 4. 

(Here P(k) is the matter power spectrum which by def- 
inition is proportional to 5'^.) In Fig[5l we plot the ra- 
tio P{k)QCDM I P{k)hcDM as a function of length scale 
of perturbation Ap. This figure implies that greater the 
length scale Ap and greater is the value of the parame- 
ter A, more is the percentage of suppression. In fact at 
Ap = lif'Mpc, matter power spectrum is suppressed by 
about 15% if the potential parameter A = 1. This im- 
plies that if A = 1 then at Ap = lO^Mpc, P(k)\cDM — 
P{k)QCDM = 0.15 X P{k)ACDM- Since P{k)ACDM at 
lO'^Mpc is at least an order of magnitude greater that 
the corresponding value at lO^Mpcfisj. it follows that 
although the percentage of suppression is larger at larger 
scales, the actual difference P{k)\cDM — P{k)QCDM is 
smaller at larger scales. 

In Fig. [6l we have plotted the perturbation in dark 
energy (5^ as a function of scale factor for A = 1. We can 
see that initially, at around z — 1000, perturbation 5rf, 
is almost zero. Once the dark energy dominated epoch 
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FIG. 9: The left plot shows the ratio of perturbation in dark energy to the matter perturbation at the present epoch for length 
scale of perturbation of 500Mpc as a function of potential parameter A. At this scale the ratio remains at almost 10""". The 
plot on the right shows the ratio at the present epoch for length scale of perturbation of 10^ Mpc as a function of parameter 
A. At this length scales, for different value of the parameter A between 0.1 and 2, the ratio varies from 0.01—0.33. 



begins the perturbations in dark energy grow. 

A. Dependence of the ratio 5^/5™ on length scales 

In Figs, m and [H] we have scaled Sm and 6^ with re- 
spect to the initial value of the metric perturbation. 
From FigU] and [6l we find that at the present epoch [z 
= 0], the ratio of dark energy perturbations to matter 
perturbations at Xp — lO^Mpc and A = 1 is given by 




Given a fixed A, the ratio 5^/5„i would depend on the 
length scale of perturbation Ap. In Fig. [3 we plot the 
value of this ratio at the present epoch as a function 
of length scale Ap. This figure shows that at small scales 
[Ap < lOOOMpc], the perturbations in dark energy can be 
neglected in comparison with the perturbations in mat- 
ter. This is because in these length scales 5^ ~ lQ~^5m- 
And since 5m itself is small, this value of 5^ corresponds 
to higher order term. In the linear regime, for scales 
Ap < lOOOMpc we can neglect the perturbations in dark 
energy (for this model) and we can treat dark energy to 
be homogeneous. The effect on the matter perturbation 
by dark energy on these scales would be through back- 
ground a{t). 

On large scales [for Ap > lOOOMpc], the dark energy 
perturbations can become comparable to 5m- In Fig. 
13 for Ap = lO^Mpc, we find that (^^/^m),,^^, = 0.17. 
Even on large scales, the perturbations in dark energy 
can be neglected if the equation of state parameter at the 
present epoch is very close to —1. This matches with the 
fact that in a pure cosmological constant model of dark 
energy with w = —1, its energy density is distributed 
homogeneously at all length scales. 

In FiglSl we have plotted the variation of the ratio 
5^ /5m at the present epoch as a function of Wo which 
is the equation of state parameter at the present epoch. 



Each value of A would result in a specific value of the 
equation of state parameter Wo as determined by the 
background equations ^ and ([5]) for a fixed set of initial 
conditions. This figure shows that this ratio 5ci,/5m ^ 
when Wo — > —1. This result is true for all length scales 
of perturbations and it is consistent with the argument 
presented in Sec|II] that perturbation in matter implies 
perturbation in dark energy if Wde 7^ — 1- 

The two plots in Fig. [51 show the ratio (5^ / 5m as func- 
tion of A at length scale of perturbations of 500Mpc and 
lO^Mpc respectively. Note that the figure on the left is 
scaled by a factor 10^. These figures imply that on large 
scales the dependence on the parameter A is stronger than 
on small scales. 



B. The role of quintessence on the matter power 
spectrum 

It is evident from Ea. (|¥T|l and Fig[5]that matter per- 
turbation is suppressed relative to ACDM model. On 
length scales Ap < lOOOMpc, the perturbation in dark en- 
ergy is negligibly small compared to the perturbation in 
matter (see Fig (7]). However, even on these scales matter 
power spectrum is suppressed relative to that in ACDM 
model. This is evident from FigO It is therefore neces- 
sary to distinguish the role of background evolution a{t) 
and perturbation in dark energy on the suppression of 
matter power spectrum relative to ACDM. 

In order to address this issue we evaluate the suppres- 
sion of matter power spectrum relative to ACDM if we 
treat dark energy as homogeneous. In such a scenario, 
this suppression would be solely a consequence of differ- 
ent background evolution relative to ACDM. Since we are 
considering quintessence model of dark energy, by ^^ho- 
mogeneous dark energy''^ we mean that we are assuming 
that quintessence field is homogeneous in the longitudi- 
nal gauge. There exists a gauge known as uniform field 
gauge where by definition the quintessence field is homo- 
geneous. However, the coordinate transformation from 
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FIG. 10: In this plot the bold line (labeled as QCDM) shows 
how matter perturbation Sm evolves with scale factor if we 
include perturbations in the quintessence field. For com- 
parison the dashed line (labeled as HQCDM) corresponds to 
Sm{a) assuming homogeneous quintessence field. In this plot 
fc = 5 X 10~''Mpc~^ and parameter A = 1. 



FIG. 11: This plot shows the ratio P{k)QCDM / P{k)HQCDM 
at the present epoch for different length scales and for two 
different value of the parameter A. In the x-axis, "Log" refers 
to the logarithm to base 10. 



uniform field gauge to the longitudinal gauge would nec- 
essarily result in non zero fluctuation in the quintessence 
field. Here our aim is to calculate the suppression of 
matter power spectrum relative to ACDM if we force- 
fully impose the assumption that the quintessence field 
is homogeneous at all length scales and compare the same 
without imposing this assumption. 

If the quintessence field is homogeneous, then the evo- 
lution of the metric perturbation $ is determined by the 
following equation : 



$ + 4-$ 

a 



3^ , ,\ a 
2 - -n„, (a) - 
a \ 2 J a 



$ = (43) 



This equation follow from Eq. (|22)) . Consequently the 
evolution of the matter perturbation is determined by 
the following equation : 



{HQCDM) ^2 ^ _ 



a ■ 
3-$ 



+ 



(a) 



}(44) 



Here the superscript "HQCDM" stands for homogeneous 
quintessence -1- cold dark matter. The above equation 
[EqlH] follows from Eqs.dH]) and (j^ . 

In Figlini we have plotted the evolution of 
{HQCDM) g^f^^-^ with the scale factor a{t). For compar- 
ison bold line in the same figure shows 
which corresponds to matter perturbations if we include 
quintessence fluctuations in our calculations. From this 
flgurc, it follows that although matter perturbation is 
suppressed relative ACDM, it is in fact enhanced in 
comparison to matter perturbation obtained by treating 
quintessence field as homogeneous. Hence perturbation 
in dark energy actually enhances matter perturbation. 
This is also intuitively expected as both dark matter and 



dark energy clusters not anti cluster in the model consid- 
ered in this paper. This is evident from FigH] and Fig [5] 
since both Sm and have same sign. The consequent 
effect of this on the gravitational potential results in this 
enhancement of matter perturbation relative to the ho- 
mogeneous dark energy case. The accelerated expansion 
of the of universe determined by a{t) dilutes the metric 
perturbation $ and consequently the overall growth of 
matter perturbation Sm is suppressed. 

In Fig. [m we plot the ratio P{k)QCDM / P{k)HQCDM 
at the present epoch as function of length scale of per- 
turbation for values of the parameter A = 0.1 and 
A = 1. On length scales Ap < lOOOMpc, we find that 
P{I^)qcdm ~ P{k)HQCDM and this is independent of 
the choice of parameter A. This implies that on these 
scales, including or excluding quintessence fluctuation in 
the perturbation equation does not influence the matter 
power spectrum significantly. The suppression of mat- 
ter power spectrum on these scales (as shown in FiglS]) is 
therefore primarily due to different background evolution 
relative to that in ACDM model. 

However, on large scales Ap > lOOOMpc, P{k)QCDM 
deviates significantly from P{k)HQCDM for larger value 
of the parameter A (see Fig fTlT) . Hence on these scales, in- 
cluding or excluding quintessence fluctuation in the per- 
turbation equation does influence the matter power spec- 
trum significantly. 

FiglHimplies that P{k)QCDM < P{k)ACDM- However, 
Figdl] impHes that P{k)QCDM > P{k)HQCDM- This 
means that P{k)HQCDM < P{k)QCDM < P{k)ACDM- 
This implies that matter power spectrum is suppressed 
relative to that in ACDM model even if treat quintessence 
field as homogeneous. This also implies that al- 
though large scale matter perturbation is suppressed 
in generic quintessence dark energy model compared 
to that in ACDM, perturbations in dark energy (in 
quintessence) enhance matter perturbation relative to 
the corresponding matter perturbation obtained by treat- 
ing quintessence field as homogeneous. This enhance- 
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FIG. 12: In the figure we plot the ratio ol dark energy per- 
turbations to dark matter perturbations for V{(fi) = l/2m^4>^. 
The results agree with those shown in Fig[7l In both the axis, 
"Log" refers to the logarithm to base 10. 



ment is significant on large scales i.e for Ap > lOOOMpc 
(see Fig. \T^. 

We compare our results with a different scalar field po- 
tential V{(jji) = im^i^^. In Fig. [T^lwe plot the ratio of 
dark energy perturbations to matter perturbations as a 
function of length scale. The parameter m is fixed to 
m — O.QAHq in natural units and the present day equa- 
tion of state parameter is w = —0.87. The results are 
consistent with those shown in Fig[7l Hence our result 
that quintessence dark energy can treated as homoge- 
neous at scales Ap < lOOOMpc is a generic result. 



perturbation in matter 5m at least in the linear regime. 
We have demonstrated this using an exponential poten- 
tial for the quintessence field. This result agrees with 
those presented in Ref.[3l[ on sub-Hubble scales. 

We have further demonstrated that quintessence dark 
energy results in suppression of matter power spectrum 
relative to KCDM model. We found that at length scale 
of Ap = lOOOMpc and for the value of the parameter 
A = 1, the matter power spectrum is suppressed by about 
4% compared to its value in the ACDM model for the 
same set of initial condition. However, at Ap = lO^Mpc, 
matter power spectrum is suppressed by about 15% com- 
pared to its value in the ACDM model. We have demon- 
strated that on scales Ap < lOOOMpc this suppression is 
primarily due to different background evolution relative 
to ACDM model. The resultant matter power spectrum 
is nearly invariant even if we assume that quintessence 
field is homogeneous on these scales. However, on much 
larger scale Ap > lOOOMpc, including or excluding fluc- 
tuation in the quintessence field results in significant 
changes in the matter power spectrum. 

All these results emphasize that dark energy can in- 
deed be treated as nearly homogeneous on scales Ap < 
lOOOMpc. However, on much larger scales (Ap > 
lOOOMpc), if the equation of state parameter deviates 
from -1, then perturbations in dark energy does influ- 
ences matter power spectrum significantly. If a definitive 
detection of perturbations in dark energy is made, it will 
certainly rule out the cosmological constant at least as 
the sole candidate of dark energy. 



IV. SUMMARY AND CONCLUSIONS 

In this paper, we have investigated the perturbations 
in dark energy. This is motivated by the fact that the 
assumption that the distribution of dark energy (with 
w 7^ —1) is homogeneous at all length scales is incon- 
sistent with the observational fact that dark matter is 
distributed inhomogeneously. On length scales compara- 
ble to or greater than the Hubble radius (Ap > lOOOMpc), 
the perturbation in dark energy can become comparable 
to perturbation in matter if Wde 7^ ~1- The model pa- 
rameters we have chosen correspond to w ~ —0.8 and 
w ~ —0.9, which are within the range allowed by Su- 
pernova observations and WMAP5 observations. Given 
this range, the evolution of perturbations differs signif- 
icantly. For scales Ap < lOOOMpc, the perturbation in 
dark energy 5^ can be neglected in comparison with the 
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